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Abstract. The Newtonian limit of fourth-order gravity is worked out discussing 
its viability with respect to the standard results of General Relativity. We 
investigate the limit in the metric approach which, with respect to the Palatini 
formulation, has been much less studied in the recent literature, due to the higher- 
order of the field equations. In addition, we refrain from exploiting the formal 
equivalence of higher-order theories considering the analogy with specific scalar- 
tensor theories, i.e. we work in the so-called Jordan frame in order to avoid 
possible misleading interpretations of the results. Explicit solutions are provided 
for several different types of Lagrangians containing powers of the Ricci scalar 
as well as combinations of the other curvature invariants. In particular, we 
develop the Green function method for fourth-order theories in order to find out 
solutions. Finally, the consistency of the results with respect to General Relativity 
is discussed. 
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1. Introduction 

The study of possible modifications of Einstein's theory of gravity has a long history 
which reaches back to the early 1920s [HEIIS]. 

Corrections to the gravitational Lagrangian, leading to higher-order field 
equations, were already studied by several authors [H [5l H] shortly after General 
Relativity was proposed. Developments in the 1960s and 1970s [71 [8l [9l [TOl [TT]. 
partly motivated by the quantization schemes proposed at that time, made clear that 
theories containing only a term in the Lagrangian were not viable with respect to 
their weak field behavior. Buchdahl, in 1962 [7j, rejected pure theories because of 
the non-existence of asymptotically flat solutions. 

The early proposed amendments of Einstein's theory were aimed at a unification 
of gravity with other branches of physics, like Electromagnetism; recently the interest 
in such modifications comes also from cosmology. For a comprehensive review, see 
[12] . In order to explain observational data, additional ad- hoc concepts, like dark 
energy/matter, are introduced within Einstein's theory. On the other hand, the 
emergence of such stopgap measures in a cosmological context could be interpreted 
as a first signal of the breakdown of General Relativity on these scales [13l [14] , and 
led to the proposal of many alternative modifications of the underlying gravity theory 
(see [T5l[16] for reviews). 

While it is very natural to extend Einstein's gravity to theories with additional 
geometric degrees of freedom, see for example [T7J [TSJ [12] for some general surveys on 
this subject as well as [20j for a list of works in a cosmological context, recent attempts 
focused on the old idea of modifying the gravitational Lagrangian in a purely metric 
framework, leading to higher-order field equations. Due to the increased complexity 
of the field equations in this framework, the main body of works dealt with some 
formally equivalent theories, in which a reduction of the order of the field equations 
was achieved by considering the metric and the connection as independent objects [21j . 
In addition, many authors exploited the formal relationship to scalar-tensor theories 
to make some statements about the weak field regime, which was already worked out 
for scalar-tensor theories [22j . 

In this paper, we shall study the Newtonian limit of fourth-order gravity theories 
in which extensions of the Hilbert-Einstein Lagrangian are considered. We are going 
to focus on the weak field limit within the metric approach. At this point it is useful 
to remind that it was already shown in 23J that different variational procedures do 
not lead to equivalent results in the case of quadratic order Lagrangians. 

In addition, we carry out our analysis in the so called the Jordan frame, i.e. 
we do not reduce the theory under consideration to a simpler one by means of a 
conformal transformation. This is due to the fact that it was already shown earlier 
[24l[25] that non- linear theories of the kind considered here, with the exception of some 
special cases, could be not physically equivalent if they have undergone a conformal 
transformation. The debate on this topic is open as can be seen in recent literature 
(see for example [li [571 [Ml [Ml [301 [31] ) ■ 

By considering some admissible choices for the gravitational Lagrangian with 
quadratic corrections, we explicitly work out the weak field limit. Such considerations 
are developed also in relation to the companion papers [32l [33l [34] where we have 
considered the Newtonian limit explicitly for f{R) gravity, taking into account generic 
analytic functions of the Ricci scalar i?, and the spherically symmetric solutions vs. 
the weak field limit respectively. 
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In principle, any alternative or extended theory of gravity should allow to recover 
positive results of General Relativity, for example in a weak limit regime, then starting 
from the Hilbert-Einstein Lagrangian 

/:o = i?, (1) 

the following terms 

Ci = R\ (2) 

C2 = Ro^pR^^ (3) 

C3 = Ro^Pf^.R"^^", (4) 

and combinations of them, represent the obvious minimal choices for an extended 
gravity theory with respect to General Relativity. Since the variational derivative of 
£3 can be linearly expressed [31 |37l [33 |3S] via the variational derivatives of Ci and 
£2, one may omit £3 in the final Lagrangian of a fourth-order theory without loss 
of generality. In this paper, we will consider the Newtonian limit of the combined 
Lagrangian ([I])-©-©, see Sec[21 as a straightforward generalization of the Einstein 
theory which is obviously recovered in low curvature regimes. 

As we said, several works focused on the cosmological implications of additional 
terms to the Hilbert - Einstein Lagrangian [38, |3S1 [HJ [32] • In particular, terms 
of the form Rr"^ (with positive n) have been taken into account to explain the 
observed accelerated behavior of the Hubble flow. Although such a term may 
lead to an alternative explanation of the acceleration effect, its singular behavior 
clearly leads to problems in the low curvature regime. The explicit exclusion of flat 
solutions is in contradiction with the basic assumptions of most of the weak field 
approximation schemes of General Relativity, the most prominent examples being the 
post-Minkowskian and the post-Newtonian approximation [32] • In addition the non- 
validity of flat solutionfQ leads to the paradoxical situation that standard linearization 
procedures can no longer be based on Minkowskian space in the lowest order and are 
therefore not well defined in such a framework. While in a purely cosmological context 
this drawback could be circumvented by basing the analysis on a curved background, 
in the lowest order, it is not obvious how such a theory could make sense at local 
scale^. In addition, it is not clear how one could have a smooth and well-defined 
transition to General Relativity for this kind of theories. Due to these drawbacks, we 
will consider here theories which allow to recover the flat solution. 

The plan of the paper is the following. Sec[2] is devoted to give the general form 
of the fourth-order field equations and their approximations at the lowest order in 
the weak field limit, i.e. the Newtonian one. Sec[3]is devoted to the solutions of the 
field equations in the Newtonian limit. In Sec[31 we develop in details the method of 
Green functions for system with spherical symmetry, while in Sec[5l we exhibit explicit 
solutions derived using the Green function. We draw conclusion in Sec[6] and present 
a possible outlook for future developments. In [Appendix A[ we present an alternative 
approach to solve the field equations, while we summarize the conventions and the 
dimensions of quantities used throughout the text in [Appendix B[ 

I As explicitly demanded by some authors, see p. 3 in 1431 . 
§ We use the term local here for the rather broad range ~ 

10^2 _ iQllm, hence encompassmg 

laboratory as well as solar-system experiments. 
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2. Gravity with quadratic Lagrangians: The field equations and the 
Newtonian limit 

In this Section, we discuss the fourth-order field equations and their Newtonian hniit. 
This higher-order, with respect to the standard second-order of Einstein field equations 
is due, as well known, to the integration of the boundary terms. These terms disappear 
in General Relativity, thanks to the Divergence Theorem, but this is not possible for 
several alternative theories of gravity, as the higher-order ones, and then the derivative 
order of field equations results augmented. 

In this paper, we are interested to achieve the correct Newtonian limit of gravity 
theories with quadratic Lagrangians in the curvature invariants. This result can be 
achieved under two main hypotheses: i) asking for low velocities with respect to the 
light speed and ii) asking for week fields. By these requests, the metric tensor is 
independent of time and second order perturbation terms can be discarded in the 
field equations (see also [32] for details). It is worth stressing that the Newtonian 
limit of any relativistic theory of gravity is related to such hypotheses and it is a 
misunderstanding to consider only the recovering of the Newtonian potential. In 
other words, a more general theory of gravity gives rise, in the Newtonian limit, to 
gravitational potentials which can be very different from the standard Newtonian one. 

In the third and in the fourth parts of this Section, we shall discuss the field 
equations in the Newtonian limit to show that the fourth-order contributions to the 
potential cannot be trivially discarded. 



2.1. General form of the field equations 

Let us now come back to the choices displayed in (II|)-([2|)-(l3])-(l4]), for which the left- 
hand side of the field equations takes the general form 

^H^^ — R^^ - -g^_iiyR, (5) 
^H^^ = 2RR^„ — -gf_i^R^ — 2R.f^^ + 2g^i,\3R, (6) 

+ 5M.^"^a/3, (7) 
'^HiJ.u — '2'R^j.al3iRu°'^'^ — ■^gfiiyRap-fSR"'^^^ + ^R^"' \(af3)- (8) 

All of the three expressions in ©-([T])-® involve fourth order differential operators. 
Due to the identity 

iff^, -4 2iJ^,+ = 0, (9) 

which holds in a four-dimensional spacetime [37], only two of the expressions in ©-([T])- 
([8]) are independent, and we are free to use any two independent linear combinations 
in our analysis. This identity gives rise to the well known Gauss - Bonnet topological 
invariant which recently acquired a lot of importance in cosmology as a possible source 
of dark energy [40j . Furthermore, for a Lagrangian comprising a general function of 
the Ricci scalar, we have 

""'"»-^«--5W-(S)^^/*.Q^. (10) 
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Here we denoted the covariant derivatives by a semicolor[j| and the d'Alembert operator 
by 110. With these considerations in mind, let us to consider the Newtonian hmit of 
such a theory of gravity. 

2.2. The Newtonian limit 

Here we are not interested in entering the theoretical discussion on how to formulate 
a mathematically well sound Newtonian limit of general relativistic field theories, for 
this we point the interested reader to[lllll51llSlll7lllHllinilSn]- In this section, 
we provide the explicit form of the field equations for the different admissible choices 
of Lagrangians collected in Introduction at the lowest, i.e. Newtonian, order. In the 
language of the post-Newtonian approximation, we are going to consider the field 
equations up to the order 0(0^^), where c denotes the speed of light. For the 
verification of our calculations we made use of a modified version of the Procrustes 
package [51]. 

We only mention, in passing, that there has also been a discussion of a somewhat 
alternative way to define the Newtonian limit in higher-order theories in the recent 
literature, see for example [52;. In this work, the Newtonian limit is identified with the 
maximal symmetric solution, which is not necessarily Minkowski spacetime in f{R) 
theories which could be singular. 

Let us start from a flat background and work out the corresponding field equations 
and hydrodynamic equations to the Newtonian order. Our conventions are that 

with a,(3 = 0,1,2,3, can be transformed to rj^p = diag(l, — 1, — 1, — 1) along 
a given curve. Latin indices i,j run from 1,2,3, the coordinates are labelled by 
a:" = (x", x^, a;^, x^) — (ci, x^, x^, a;'^) . We start with the following ansatz for the 
metric 

ffoo = l-^ + 0(c-^), 

90a = 4r/lOa + C(c"^), 

gab = - (1 + 5at + 0{C-^). (11) 

Apparently, the orders involved in this ansatz for the line element reach beyond the 
Newtonian order, which we are mainly interested in this work. 

On the matter side, i.e. right-hand side of the field equations, we start with the 
general definition of the energy-momentum tensor of a perfect fluid 

Tal3 = {PC^ +Ilp+p) UaUp - pgaf3, (12) 

here H denotes the internal energy density, p the energy density, and p the pressure. 
Following the procedure outlined in [S3], we derive the explicit form of the energy- 
momentum as follows 

Too = + O (c-2) , (13) 

Toa = cpv'^ + O (c-i) , (14) 

Tab ^ pv^'v'' + pSab + O {C-^) . (15) 
II We denote partial derivatives with respect to the coordinates by a comma. 
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T 



(16) 

with a generalized tensor H^i,, being a combination of the expression^ specified in 
©-([B])-©-®, which, in turn, depend on the final form of the Lagrangian. 



2. 3. The quadratic Lagrangians and the Newtonian limit of the field equations 

Let us consider now the field equations, in the Newtonian limit, for the possible 
quadratic Lagrangians which we compare to the Newtonian limit of the standard 
Hilbert - Einstein Lagrangian. It is important to stress that the field equations, in 
Newtonian limit, are considered up to the order 0(c~^) while the vector component 
to the order 0{c~^) is related to the post-Newtonian limit of the theory (see [S3] 
for details). Up to the Newtonian order the left-hand side of the field equations, i.e. 
dS])-©-©-®, takes the following forn^ for the metric given in pTjl : 



The Hilbert - Einstein Lagrangian (£o — R) 



Hab — 



The R"^ - Lagrangian {Ci ~ R? 



li/oo = ^ (2 W - V^C/) 

^a^a — 0, 



{V^U - 2VV) 5ab 
The RapR"^ - Lagrangian (£2 = RapR"^^) 



Hoa = 0, 



1 r 



{V^U - SV^V) Sab + (3V2y - v^c/) 



The RafjjsR"'^^'^ - Lagrangian (£3 — Ra/s-ysR 



Hoa — 0, 
'Hab - - 



Obviously "-ffpi^ is the Einstein tensor. 

Here we made use of the following operator definition 



(17) 
(18) 
(19) 

(20) 
(21) 
(22) 

(23) 
(24) 
(25) 

(26) 
(27) 
(28) 



as well as V" := V^V^ 



The Newtonian limit of metric gravity theories with quadratic Lagrangians 



7 



2-4-. The combined Lagrangian 

Let us now combine the different terms of the last section into the same Lagrangian. 
This combination is the basis for our investigation for the remaining part of the paper. 
Since terms resulting from i?" with n > 3 do not contribute at the order O (c~^) , the 
most general choice for the Lagrangian is 

C = a^R + a2R^ + b^R^.R^" + c^R^p^,R"P^'r (29) 

Due to the identity given in it is sufRcient to study 

C = aiR + a2R^ +biRf,^R^"', (30) 

in four dimension, where we introduced the constants ai, 02, foJJJ- If we take into 
account the results from Unil-lIlSl) as well as (IT7l)-(l25l) the explicit form of the field 
equations (fT6|) up to the Newtonian order is 



- 2aiVV + (802 + 2hi)V'^V - (4a2 + 2h{)V^U = SnGp , 



(31) 



ai(VV - V^C/) - (802 + 361)^^ + (4a2 + bi)\/'^U 
+ (8a2 + 3bi)\7^V - (4a2 + bi)V^U + ai{U ~ V) 
the equations which we are going to solve. 



Sab 

= 



(32) 



3. Considerations on the field equations in the Newtonian limit 

In this Section, we are going to formulate the problem to solve the field equations 
(|3ip - (|32|) in the most general way. It is worth noticing that the isotropic coordinates 
for the metric pip allow to search for solutions independently of the symmetry of 
the physical system (which can be spherical, cylindrical etc.). The results which 
we are going to achieve are completely general since we will search for solutions in 
terms of Green functions. However, being the combined Lagrangian pop built up 
by various terms, the field equations strictly depend on the coupling constants. As 
we will see, the value of such coefficients have a crucial role for the validity of the 
approach since, from a physical viewpoint, we have to obtain the Newtonian limit of 
General Relativity as soon as the quadratic corrections disappear. This aspect of the 
problem is not accurately faced in the literature and can lead to wrong conclusions. 
Our aim is to develop a method which allows to control, step by step, the Newtonian 
limit in agreement with the results of General Relativity. This is possible at three 
levels: Lagrangian, field equations and solutions. In the second part of this Section, 
we will analyze the various cases of field equations considering particular values of 
the coefficients. Specifically, we will take into account the values where the proposed 
approach fails. It is interesting to note that any time the Hilbert - Einstein term 
is absent into the Lagrangian, the field equations are fourth order (see Table [Ij. In 
other words, if the Hilbert - Einstein term is not present, we do not recover the 
Laplace/Poisson equations. 

tt Note that [ai] = [length]", [02] = [h] = [length]^. 

tt The coefficients of ^ are different from ones of Il30|l . 
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By introducing two new auxiliary functions {A and -B), the equations (|5T|) - (P^ become 



^2) ^02+^1 A 

2a2 + bi 
2a2 



ai 



B 



2a2 + bi 
= SnGp , 



26i(3a2 + 
ai(2a2 + bi) 



A 



2a2 + hi ' 

where A and B are Hnked to U and V via 

A:^a^{V^U), 

B := (4a2 + bi)V - {802 + 3bi)U . 



(33) 
(34) 

(35) 
(36) 



Obviously we must require ai(2a2 + 61) ^ 0, which is the determinant of the 
transformations (P5)) - (|36p . Let us introduce the new function <i> defined as follows: 



A + V^B. 



(37) 



At this point, we can use the new function $ to decouple the system In fact 

we obtain 



26i(3a2 + 6i) 



where r/ 



4a2 



ai(2a2 + hi) 
+ 5i , 26i(3a2 + 6i) 



602 + bi 
2a2 + bi 



V^B 



2a2 + bi 



V^B 



2a2 + 61 * ' ai(2a2 + 61) 
in terms of the Green function 5/(x, x') defined by 



(38) 
(39) 

V^$. We are interested in the solution of (l38t 



where 



B(x)=y, J dVg/(x,x')a/(x') 



a7(x) :=87rGp(x) - VV7(x), 



(40) 



(41) 



and Yj being a constant, which we introduced for dimensional reasons. Then the set 
of equations ([HTjl - iP^ is equivalent to 



26i(3a2 + 6i) e , 602 + &i 4 , 

— - — -— V°a/(x,x ) + - — --_V^e/(x,x ) 
ai(2a2 + Oi) 2a2 + Oi 

- ^^v^c?,(x, x') = -r,-i5(x - x') , 

2a2 + oi 
V2$(x),5,b-$(x),,b = 0, 



(42) 
(43) 



where (5(x — x') is the 3-dimensional Dirac (5-function. The general solutions of 
equations ((5T1) - ([5^ for J7(x) and V{x.), in terms of the Green function Qj{x,x') and 
the function $(x), are 
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87rGp(x') 



2a2 + bi 
8a2 + 35 



ai(2a2 + 61) 



2ai(2a2 + 5i) 



(44) 



v(x) = 



(4a2 + bi)V^ - ai 
bi) 



2ai(2a2 
4a2 + 5 
2a2 + bi 
4a2 + 61 



iv2,$(x')- 



(i^x'CJ/(x, x') 

2bi(3a2 + 61) 
ai(2a2 + 5i) 



87rG'/9(x') 
V^,ci>(x') 



r*(x) 



2ai(2a2 + 6i)^'"'' ^"^^^ 

Equations PT|) - ([5^ represent a coupled set of fourth order differential equations. 
The total number of integration constants is eight. With the substitution ((37)) . it 
has been possible to decouple the set of equations, but now the differential order 
is changed. The total differential order is the same, indeed we have one equation 
of sixth order (|38p . and another equation of second order (|39l) . while previously we 
had two equations of fourth order. Obviously, the number of integration constants is 
conserved. The possibility to decouple the field equations ((3T|) - ([32]) is strictly realted 
to the choice to express the auxiliary field A in terms of B by inverting the relation 
([57]) . deriving Eq. (|5T|) and reducing to the second order. In [Appendix A[ we 
discuss a different where the set of equations ([5T|) - ([5^ remains of forth order by 
using the relation between B and A. 



3. 2. Field equations for particular values of the coupling constants 

In this subsection, we want to analyze the behavior of the field equations ((3T|) - ((32)) 
for those values of the coupling constants oi, 02, 61 where the transformations (|35p - 
(|36p do not hold. Specifically, in Table (TJ we display several cases of the field equations 
pip - ([5^ . for different choices of the coupling constants, where the determinant of 
transformations - is zero. First of all, we have to note that, for 02 = &i = 
(Case i in the Table [1]), we trivially obtain the same result of General Relativity in 
isotropic coordinates. Furthermore, by asking for U = V, the spatial equation is 
satisfied. It is straightforward to derive as solution the Newton potential (see § 15.31 
for details). Particularly interesting is also Case iv, where both terms Ra/sR'^^ and 
R^ give similar contributions. In other words, we have the same situation of the 
Lagrangian £ = aiR + a2Ra0R°'^ with a redefinition of the couplings. However, this 
Lagrangian is compatible with the transformations ((35|) - ((36)) and then the results of § 
l3.1l hold. In Cases (ii - iii - v - vi - vii), the differential operator never appears as a 
linear term since the invariants R^ and RapR"'^ give rise to higher order terms in the 
field equations. In these cases, the full field equations (not in the weak field regime) 
give g^i,n\R — R-^^y for the Lagrangian R?' and — 2i?(^|".|t/)Q + ni?pj/ + g^y 

R''^;c.p for 

RapR"^ which are fourth order equations. In the weak field regime, one obtains the 
equations reported in Table [T] As we will see in § 15.31 devoted to the solutions, all 
these cases do not present a Newtonian potential. 
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Cases 


Choices of oi, 02, 61 


Corrcsponding field equations 


i 


a2 = 
61 = 


v-u 




V ~U =0 


ii 


oi = 
61 = 


V4( 


2y - [/) = ^ 
V2(2y- [/) 


Sab - 


V2(2V^- [/) 


= 


iii 


01 = 

02 = 


V'*( 


c/-y) = -^ 

V2(t/-3V^) 


fci f^' 

Sab - 


V^{U -W) 


= 

.ab 


iv 


bi = -2a2 


V2 


V^y-fliV^y = 47rGp, 
ai(y-[/)-2a2V2(F-[/) 

ai(V^-t/)-2a2V2(V^-f7) 


(Job 

= 

ab 


V 


oi = 

bi = -4a2 






Sab - 




= 


vi 


ai = 
61 = -2a2 






Sab ~ 


W^{V -U) 


= 0; 


vii 


oi = 
bi - — 3- 


V2 


2V + 


U) = ^ 




= 

,a6 



Table 1. Explicit form of tiie field equations for different choices of the coupling 
constants for which the determinant of the transformations l|35| l- l|36|l vanishes. 
Cases i, ii, iii are the Lagrangians introduced in the Sec l2.3l (R. , i?,^^/?'"^). 



4. Green's functions for spherically symmetric systems 

We are interested in the solutions of field equations at order 0{c~^) by using the 
method of Green functions. We have to stress that the method of Green's functions 
does not work in the general case since the field equations are non-linear. However, 
the Newtonian limit of the theory (based on the hypothesis that metric nonlinear 
terms can be discarded) allows that also the field equations result linearized). By 
solving the field equations with the Green function method, one obtains, as a first 
result, the solution in terms of gravitational potential in the point-mass case. Then by 
using the equations (|44l) - (|45)l . obtained in the weak field limit and then in Newtonian 
linear approximation for a spatial distribution of matter, we obtain, in principle, the 
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gravitational potential for a given density profile. If the matter possesses a spherical 
symmetry, also the Green function has to be spherically symmetric. In this case, the 
correlation between two points has to be a function of the radial coordinate only, that 
is: Q{x,k') = Q{\x — x'|). It is important to stress again the fact that the approach 
works if and only if we are in the linear approximation, i.e. in the Newtonian limit. 

4.I. A general Green function for the decoupled field equations 

Let us introduce the radial coordinate r := |x — x'|; with this choice, equation 
for r ^ becomes 

26i(3a2 + bi)Vlgiir) + ai(6a2 + 6i)V^5/(r) - ajVlGiir) = , (46) 

where = r^^dr{r^^dr) is the radial component of the Laplacian in polar 
coordinates. The solution of (l46l) is: 



Giir) = Ki,i - - 
r 




where Kj i, Ki^2, Kj ^, -ft'/, 4, i^^/.s, Kj q are constants. The integration constants 
Ki i have to be fixed by imposing the boundary conditions at infinity and in the 
origin. A physically acceptable solution has to satisfy the condition Q{x,x') — > if 
|x — x'l — > oo, then the constants Kj i, Ki^4, Ki q in equation ([17|) have to vanish. 
We note that, if 02 = 6i = 0, the Green function of the Newtonian mechanics is 
found. In this case, we have the complete analogy with the Electromagnetism. More 
precisely, when we do not consider higher-order tems than Hilbert - Einstein one in the 
gravitational Lagrangian, we obtain, in the Newtonian limit, a field equation analog 
to the electromagnetic one for the scalar component (electric potential). This means 
that we have the same form of the Green function [55] . 

To obtain the conditions on the constants Kl^2^ Kj^^ we consider the Fourier 
transform of Q{x,x.'): 

^^C7z(k)e^'^-(— ). (48) 

^//(x, x') depends on the nature of the poles of |k| and on the values of the arbitrary 
constants ai, 02, bi. If we define two new quantities Ai, A2 G TZ: 



Oi 2.{Aa2 + Oi) 



we obtain: 



e/(x,x') = 




2 |x— X 



(50) 



This Green function corresponds to the one in (|T7l) . Obviously, we have three 
possibilities for the parameters Ai and A2. In fact, Ai and A2 are related to the 
algebraic signs of ai, 02, hi and then we can also achieve real values for such 
parameters. This means that we have three possibilities: both imaginary, one real 
and one imaginary. In Table [2l we provide the complete set of Green functions 
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5/(x, x'), depending on the choices of the coefhcients 02 and bi (with a fixed sign 
of ai). The various modahties in which we obtain the Green functions are due to the 
various sign combinations of the arbitrary constants. In general, the parameters Ai^2 
indicate characteristic scale lengths where corrections to the Newtonian potential can 
be appreciated. It is worth noticing that, thanks to the forms of the Green functions 
(see Table [5]), the Newton behavior is always asymptotically recovered. When one 



Cases 



Choices of 02, 61 



Green function C?/(x, x') 



5i < 

3a2 + 61 > 



At A9 



-^21=' 



5i > 



802 + 5i < 



18 |x-x'| 



cos( Ai |x— x' I ) cos(A2 |x— x' I ) 



A, A„ 



A? 



A? 



5i < 

3a2 + bi <0 



18 IjT 



cos(A2|x— x'l) 
A? 



Table 2. The complete set of Green functions for equations II42I I. The scale 
lengths are: Ai := \ai/bi\^^'^, A2 := |ai/2(3a2 + bi)\^^^. It is possible to have a 
further choice for the scale lengths which turns out to be dependent on the two 
knows length scales. In fact, if we perform the substitution Ai ^ A2, we obtain 
a fourth choice. In addition, for a correct Newtonian component, we assumed 
ai > 0. In fact when a2 = bi = the field equations JSTJ and II32I I give us the 
Newtonian theory of gravity if ai = 1. 



considers a point-like source, p oc (5(x), and by setting $(x) — the potentials (|44 
(|45|) are proportional to C//(x,x'). Without losing generahty we have: 



C/(x) - 

X 



= -Ai |x 



D-A2 |x 



Ui- 



U2 



(51) 



where Uq, Ui, U2 are some integration constants. An analogous behavior is obtained 
for the potential V^(x). We note that in the vacuum case we found a Yukawa-like 
corrections to Newtonian mechanics but with two scale lengths related to the quadratic 
corrections in the Lagrangian pO|) (see also the above expressions (149]) ). This behavior 
is strictly linked to the sixth order of (|^^ . which depends on the coupled form of the 
system of equations (|3ip - ([5^ . In fact if we consider the Fourier transform of the 
potentials U and V: 



C/(x) 



u(k) e' 



(27r)3/2 

the solutions of equations ([3T|) - p2|) are 

S'k 47rG'[ai 



y(x) 



d^k 

(27r)3/2 



)(k) 



C/(x) 

y(x) 



(8a2-H36i)k2]p(k)f 



(27r)3/2 k2(ai - 6ik2)[ai + 2(3a2 + foi)k2] ' 

d^k AnG[ai + (4q2 + bi)k2]p(k)e''^-^ 
(2^)3/2 k2(ai - 6ik2)[ai + 2(3a2 + 6i)k2] ' 



(52) 

(53) 
(54) 
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where /5(k) is the Fourier transform of the matter density. It is possible to show, by 
applying the Fourier transform to the potentials U and V, that the poles in Eqs. ([55)1 
- ([M)) are always three. 

. . M 

Finally, if p(k) = - — — -:- (the Fourier transform of a point-like source) the 

solutions ([551) - ([51)) are similar to (|5T|) . In fact, if we suppose that 61 7^ and 
3a2 + 61 ^ 0, the solutions ^ - (HH) are 

f/(x) ^^(1- ^e-^^l-l + le-^^^A , (55) 



ai|x| \ 3 3 

4-2. Green functions for particular values of the coupling constants 

The obtained Green functions deserve some comments. First of all, we have to consider 
the particular values of the parameters where the general approach developed in § 14.11 
does not work. For example, if bi = 0, we have only one Yukawa-like correction. The 
Green function have to satisfy the equation 



3Vie/(x,x') ~ -^V2e/(x,x') = -Vi<5(x-x'), (57) 

obtained from the (|42|) by setting 61 = 0. In this case, the Green function (Fourier 
transformed), is: 

g/(k)-^ , ^^f/\. , (58) 

and the Lagrangian becomes: L — aiR + a2R^ ■ Since at the level of the Newtonian 
limit, as discussed, the powers of Ricci scalar higher then two do not contribute, we 
can conclude that ([55)1 is the Green function for any f{R) - theory at Newtonian order, 
if f{R) is some analytical function of the Ricci scalar. The same result is achieved 
considering a particular choice of the constants in the theory, e.g. &i = —2a2- In 
Table [T] (case iv) , we provide the field equations for this choice and the related Green 
function is: 

g(2a.V^-..V^)(k)0C ^^^^^/^^^^^^ . (59) 



The spatial behavior of ((58l) - ((59)) is the same but the coefficients are different since 
the theories are different. The interpretation of the result is the same than that in 13.21 
since we have to take into account a proper scale length. In fact Eq. (|55)) presents a null 
pole for k^ = which gives the standard Newtonian potential, and a pole in k^ — \\, 
which gives the Yukawa - like correction. Finally, we need the Green function for the 
differential operator V*. From Table [U the field equations present always a quadratic 
Laplacian operator (Case i excluded). This means that the equation to solve is: 

V4g(v^)(x',x) =5(x-x'). (60) 
By introducing the variable r = |x — x'| 7^ 0, we have that Eg. ijSn)) becomes 
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^tGiv-){r) = (61) 

with solution 

g(V4) (r) = Ki^i + ^ + Ki^sr + Ki^rr^ , (62) 

where Kj^s, Ki^7 are generic integration constants. 

Let us now consider the fact that the Green function has to be nuU at infinity. 
The only possible physical choice for the squared Laplacian is: 

g(V4)(x,x') oc (63) 

Considering the last possibility, we will end up with a force law increasing with distance 
[To] . In conclusion, we have shown the general approach to find solutions of the field 
equations by using the Green functions. In particular, the vacuum solutions with 
point-like source have been used to find out directly the potentials, however it remains 
the most important issue to find out solutions when we consider systems with extended 
matter distribution. 

5. Solutions by the Green functions in spherically symmetric distribution 
of matter 

In this section, we explicitly determine the gravitational potential in the inner and 
in the outer region of a spherically symmetric matter distribution. This is a delicate 
problem since the Gauss theorem is not valid for the gravity theories which we are 
considering. In fact, in the Newtonian limit of General Relativity, the equation for 
the gravitational potential, generated by a point-like source 

^iGNew.mech. (x, x') = -47r(5(x - x') (64) 

is not satisfied by the new Green functions developed above. If we consider the flux 
of force lines 'P New.mech. defined as 



New.mech. • i , i q 



= -GMV:^g New.mech. (x, x') , (65) 

1^ ~ ^ r 

we obtain, as standard, the Gauss theorem: 

dT, F New.mech. ■ n(x M , (66) 



s 

where S is a generic two-dimensional surface and n its surface normal. The flux of fleld 
^New.mech. On the surface E is proportional to the matter content M, inside to the 
surface independently of the particular shape of surface (Gauss theorem, or Newton 
theorem for the gravitational field [56]). On the other hand, if we consider the flux 
defined by the new Green function, its value is not proportional to the enclosed mass 
but depends on the particular choice of the surface: 

F New.mech. ' f^ oc Ms . (67) 

s 
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Hence is a mass- function depending on the surface E. Then we have to find the 
sohition inside/outside the matter distribution by evaluating the quantity 

ydVg,(x,x')p(x'), (68) 
and by imposing the boundary condition on the separation surface. 

5.1. The general solution by the Green function C//(x, x') 

By considering the expressions (|44p and ([IS)) with the Green function ([50]) and by 
assuming $(x) = 0, we have 

. 4.Gr.i^^m^lM^ / d3x'g,(x,x')p(x') , (69) 
ai(2a2 + 0i) J 



y(x) =4.Gr.(l^i±^lM^ / '^^x'e.(x,x')p(x') 
ai(2a2 + 6i) J 



(70) 



We have to note that the hypothesis, <i>(x) = 0, is not particular, since when 
we considered the Hilbert-Einstein Lagrangian to give the Newtonian solution, we 
imposed an analogous condition. In fact, considering the spatial components of the 
Einstein equations and then, by asking for $ = 0, we get the condition that the two 
metric potential U and V have to be equal (see Case i in Table [T|) . In the general 
case, the role of $ is given hy V — U. From the time - time component, we can obtain 
an expression for U. The next step, in principle, is to search for a non-trivial solution 
This task is very difficult in general but can be realized for some particular cases. 
For example, since the Green function can be found under the spherical symmetry 
hypothesis, also the spatial distribution of matter have to be spherically symmetric. 
Denoting the radius of the sphere with total mass M by ^, we have the matter - 
density function 



p(x) = ^e(e-|x|), (71) 



3M 

where Q{£, — |x|) is the Heaviside function. For the potential U{x.), we obtain the 
implicit expression 



X / d|x'||x'|2 /" d<j)' [ de'sme'gi{y:,x'), (72) 
Jo Jo Jo 

and an analogous relation for V(x). After some algebra, we get the explicit form 



C^»n(x) 



(27r)3/2 GM 



Aj(2 + 3Aie) -8A: 

AjA2 



^ Ix|2 



+ 8e-^^^(l-fA,0 ""^^t^'"'^ 



Af|x| 



(73) 
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oi |x| 

4(27r)3/2 GM 



[AiCcosh(AiO -sinh(AiO]- 



ai Xie' ' '''' |x| 

+ ,3,3 [A2^cosh(A2e) - sinh(A20]^^. (74) 

ai |X| 

The relations ((73 |) -([74 | give the solutions for the gravitational potential U inside and 
outside the constant spherically symmetric matter distribution. A similar relation is 
found for V{k). The boundary condition on the surface |x| = ^ is satisfied: 

UUO-UoutiO ^ 0- (75) 

We note that the corrections to the Newtonian terms are ruled by ^/(x, x'). In order 
to achieve the behavior of the potential inside the matter distribution, let us perform 
a Taylor expansion for Ai|x| ^ 1. We have: 

sinhfAlxl) , ,n , , 

~ constant + xp + ... . 76) 

A|x| 

Like for the standard Newtonian potential, this means that the inner solution of the 
corrected potential is traced by the matter distribution. 

The outer solution has to be discussed in detail. For fixed values of the distance 
|x|, the external potential Uout{^) depends on the value of the radius ^, then the Gauss 
theorem does not work also if the Bianchi identities hold [11]. In other words, since the 
Green function does not scale as the inverse distance but has an exponential behavior, 
the Gauss theorem ()67p does not hold. This means that the potential depends on the 
total mass and on the matter - distribution in the space. In particular, if the matter 
distribution takes a bigger volume, the potential \ Uout{^)\ increases and viceversa. We 
can write 

A^ cosh(AO - sinh(AO .„^ 
hm — — = oo ; (77) 

obviously the limit of ^ has to be considered up to the maximal value of |x|. The 
term defined in (ffT]) can be defined as a sort of geometric factor which takes into 
account the spatial matter-distribution. The limit puts in evidence the dependence 
on the matter distribution of the outer potential. The spherical symmetry allows to 
find out the Green functions but, in principle, they can be achieved also without this 
hypothesis. At this point, it is interesting to consider the physical meaning of the 
Green functions. Eq. ((50| represents the correlation (i.e. the interaction) between two 
points in the space, that is, it gives the possibility to calculate the potential in a given 
point as a function of the "charge" (the point - mass) of another point. Such a charge 
is described by a Dirac delta function which is the source term of the Green function. 
By summing up the contributions of all infinitesimal volume elements, we obtain the 
potential for a given matter - distribution. This analysis can be concluded with some 
considerations related to the behavior for |x| ^ ^. This means that we are moving 
away from the matter distribution. Such a limit can be given also as ^ — > 0, that is: 

^.^3Aecosh(Ag)-sinh(A0^^ (78) 
5^0 A3^3 
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For Uouti^), we have 



hm [/o„t(x) = 



5^0 ai |x| 3ai |x| 

3ai |x| 

We can choose ai — (27r)"^/'^, and then, for a point - hke mass, we have: 



(79) 



, , GM 4GAfe-^il''l lGMe-^^\''\ 

nn [/o„t X - — - ^ + . 80 

e^o |x| 3 |x| 3 |x| 

The last expression is compatible with the discussion in §. 14.11 
5.2. Further solutions by the Green functions t//(x, x') 

For the sake of completeness, let us derive the expression for the potential [/(x) for 
the other two Green functions in Table [2l By performing a similar calculation, but 
now considering Case ix in Tabled we obtain 



C^jn(x) 



GM I A|(3A|e^-2^f8A| _ ^^^^ 

sin(Ai |x|) 



--2[cos(Aie) + Aiesin(AiO] . , , 
Af Ai|x| 

+ ^[cos(A2C) + A2Csin(A2e)]^^^^^ \ , (81) 

GM 4(27r)3/2 GM ^ . 

f/o«t(x = — -3-3 sm AiO 

|x| ai AfC^ 

, , ,,,cos(Ai|x|) (27r)3/2 GAf ^ . 

- Ai^cos Aig)] + sm AzO 

|x| ai A^f'^ 

\ c /\ eMC0s(A2|x|)l 

-A2Ccos(A2C)] • (82) 

Also in this case the boundary conditions (|75p are satisfied. The considerations of 
preceding subsection hold also for the solutions (|5T|) - ([5^ . The only difference is that 
now we have oscillating behaviors instead of exponential behaviors. The correction 
term to the Newtonian potential in the external solution can be interpreted as the 
Fourier transform of the matter density p(x). In fact, we have: 

d'^' , r, -.k.x' _ 3M sin(|k|0-|k|ecos(|k|0 



(27r)3/2'^' ' (27r)2/3 |k|3^3 

and in the point - like mass limit, it is 

we obtain again the external solution for point-like source as limit of ((82|) : 
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, , GM 2GAfcos(Ai|x|) , 1 GMcos(A2|x|) , , 

lim [/o„t(x) = — — - h- . (85) 

e^o |x| 3 |x| 6 |x| 



3 |x 

Finally for the last case in Table ^ we have 



C/.„(x) 



GM f Aj(3A^g^- 2) ~8A, 
AfAl 

+ Af' + Ai|x| 



-2[cos(A2e) + A2Csin(A20] 

At 



sin(A2|x|) 
A2|x| 



(86) 



Uout (x) 



GA/ 4(27r)3/2 GM 

X [Ai^ cosh(Ai^) - sinh(Ai^)]- 



-Ai|x| 



(27r)3/2 GM 

w 



ai 



[sin(A20-A2^cos(A20] 



cos(A2|x|) 



(87) 



The limit of point-like source is valid also in this case, that is: 



lim Uouti^) 



GM 4GA/e-^il''l 
1^~3 bd 



1 GM cos(A2|x|) 
3 bd 



The results ([5Dt - ([55t - ([551) means that, for suitable distance scales, the Gauss 
theorem is recovered and the theory agrees with the standard Newtonian limit of 
General Relativity. 



5.3. Other solutions and their physical consistency 

In Table [31 we provide solutions, in terms of the Green function of the corresponding 
differential operator, for the field equations shown in Table [T] Case i corresponds to 
the Newtonian theory and the arbitrary constant ai can be absorbed in the definition 
of matter Lagrangian as above. The implicit solution is: 



[/(x) = y(x)^G/dV^(^^. 

J |x-x'| 



(89) 



For Case iv, we have: 



C/(x) = F(x) = G / d-^x 



p(x') 



(90) 



The solutions make sense only if ai/a2 > 0, which gives a scale- length. Also in this 
case, we can have different signatures for ai and 02 which give oscillating corrections to 
the Newtonian potential. We have to note that both above cases have been solved with 
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the hypothesis <i> = 0. These two cases are the only ones which exhibit the standard 
Newtonian Umit (obviously the former) . The remaining cases can exhibit divergences 
and incompatibihties. This is obvious since, as discussed in § 13.21 the absence in the 
Lagrangians of terms Unear in the Ricci curvature scalar gives field equations with 
higher - order Laplacian operators (see Cases ii, iii, v, vi, and vii). Precisely, without 
terms like V^?7 + ... = p, we could not achieve regular Newtonian-like behaviors. 
This fact could give problems in comparing inner and outer solutions with respect 
to matter distributions. In other words, the Newtonian potential is necessary not 
only to achieve physically interesting situations but also, from a mathematical point 
of view, to regularize solutions. In fact, Case ii presents an incompatibility between 
the solution obtained from the 00 - component and the one from the a6-component. 
The incompatibility can be removed if we consider, as the Green function for the 
differential operator V"*, the trivial solution: ^?(v*)|b — const. Only with this choice, 
the arbitrary integration constant Uq can be interpreted as GM. However another 
problem remains: namely the divergence in the origin and then we can conclude that 
the solution 

CM 

2V{^)-U{^) = - (91) 

|x-x'| 

holds only in vacuum. 

Besides, terms like / d'^x'C/fvi) (x, x')p(x') have to be discussed for the choice 
(|63p . The field equation with (see Table [1]) gives 

V^t/(x) cx '^'x'^^ = -47rV^p(x) ^ -47rp(x); (92) 

which is consistent only if p(x) = 0. Due to these considerations, also in the remaining 
cases, we can consistently consider only vacuum solutions. 

6. Conclusions and outlook 

In this paper, we have studied the Newtonian limit of gravitational theories whose 
action presents quadratic curvature invariants beside the standard Ricci curvature 
scalar of General Relativity. In particular, we have considered the problem to find out 
solutions of the field equations developed up to the perturbation order O (c^^). This 
is intended as the Newtonian limit while, taking into account terms up to O (c^^) and 
beyond is the post-Newtonian approximation (see for example [54]). 

After deriving the full fourth - order field equations, we have developed the metric 
and the stress - energy tensors in the Newtonian limit. The main metric quantities, in 
this limit, are the two gravitational potentials U and V which are the solutions of the 
field equations both in presence and in absence of matter. At the order c^^, quadratic 
curvature invariants give rise to and operators acting on U and V in the field 
equations. Our task has been to develop an approach to solve such equations and to 
find out corrections to the Newtonian potential emerging, as standard, from General 
Relativity. 

The method consists in searching for suitable combinations of the gravitational 
potentials U and V by which it is possible to decouple the field equations. After 
field equations are suitably decoupled, one can define Green's functions which allow 
to obtain the potentials. These potentials, however, strictly depend on the coupling 
parameters appearing in the Lagrangian of the theory. Conversely, such coupling 
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Cases 


Solutions 


Newtonian behavior 


i 


C/(x) - V{^) - ^ /d^xY'"'),, 


yes 


ii 


2V{^) - U{^) = §1 

2Vi^) - t/(x) = ^/d3x't;(v.)(x,x')p(x') 


no 


iii 


Ui^) = ^ - ^/d3x't;(v4)(x,x')p(x') 
F(x) = f|-^/dVg(v4)(x,x')p(x') 


no 


iv 


C/(x) = AttG J d3x'g(2,,v4-a,v^)(x,x')p(x') 
F(x) = 47rG / d3x'g(2a,v4-aiv^)(x,x')p(x') 


yes 


V 


= Til 


no 


vi 


W - ^ + ^ / rf'x'CJf V.) (x, x' V(x') 
F(x) = 3^/dVe(v^)(x,x')p(x') 


no 


vii 


C^(x) = i 

= -^R + ^/d^x't;(v.)(x,x')p(x') 


no 



Table 3. Here we provide the solutions of the field equations in Table [T] The 
solutions are found by setting 'I'(x) = in the ab - component of the field equation 
l|43| l or l)A.7|l . The solutions are displayed in terms of the Green functions. Uo is 
a generic integration constant. 



constants allow to classify the field equations, and then the solutions, selecting, in 
particular, some singular cases. 

A detailed discussion has been developed for systems presenting spherical 
symmetry. In this case, the role of corrections to the Newtonian potential is clearly 
evident. In general, such corrections are oscillating behaviors or Yukawa-like terms. 
This means that one of the effects to introduce quadratic curvature invariants is 
to select characteristic scale lengths which could have physical interests as we will 
discuss below. Besides, such corrections invalidate the Gauss theorem because any 
matter distribution depends on such scale lengths. Furthermore, if the Newtonian 
potential term is not present, there could be compatibility problems and some solutions 
are physically consistent only in vacuum. Furthermore, for spherically symmetric 
distributions of matter, we discussed the inner and the outer solutions and the 
boundary conditions. 

From a physical viewpoint, this systematic work is needed in order to fully 
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develop the weak field limit of such relativistic theories of gravity and then compare 
them with observations and experiments. In fact, the correct interpretation of data 
strictly depends on the self-consistency of the theory and, viceversa, data correctly 
interpreted could definitively confirm or rule out deviations from General Relativity 
[57] , It is worth pointing out that extended or alternative theories of gravity seem 
good candidates to solve several shortcomings of modern astrophysics and cosmology 
since they could address several issues of cosmological dynamics without introducing 
unknown forms of dark matter and dark energy (see e.g. [HI [38]). Nevertheless, a 
"final" alternative theory solving all the issues has not been found out up to now and 
the debate on modifying gravitational sector or adding new (dark) ingredients is still 
open. Beside this general remark related to the paradigm (extending gravity and/or 
adding new components), there is the methodological issue to "recover" the standard 
and well-tested results of General Relativity in the framework of these alternative 
schemes. The recovering of a self-consistent Newtonian limit is the test bed of any 
theory of gravity which pretends to enlarge or correct the Einstein General Relativity. 

Taking into account also the results presented in [32l [33] , it is clear that only 
General Relativity presents directly the Newtonian potential in the weak field limit 
while corrections (e.g. Yukawa-like terms) appear as soon as the theory is non-linear 
in the Ricci scalar. This occurrence could be particularly useful to solve the problem 
of missing matter in large astrophysical systems like galaxies and clusters of galaxies 
as discussed in [T4j [58] . In fact dark matter (and dark energy) could be nothing else 
but the effects that General Relativity, experimentally tested only up to Solar System 
scales, does not work at extragalactic scales and then it has to be corrected. Assuming 
this alternative point of view, we do not need to search for unknown ingredients, up 
to now not found at fundamental level, but we need only to revise the behavior of 
gravitational field at infrared scales. These scales could be ruled by corrections to the 
Newtonian potential, as shown in this paper. In forthcoming researches, we intend to 
confront such solutions with experimental data, as done in [Ml [58] , in order to see if 
large self-gravitating systems could be modelled by them. 
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Appendix A. Alternative approach to solve the field equations 

In this Appendix, we discuss an alternative approach to solve the field equations where, 
instead of using the relation A = $ — V^-B as in § 13.11 to obtain the solutions (|44|) - 
(|45p. we adopt the inverse relation V^B = ^ — A. As noted above, the first relation 
makes the differential degree of system increase but we have a relation between the 
solutions A and B. In the second case the differential degree remains the same but 
we have a non-local relation between the solutions: 



,A(x') -$(x') 



(A.l) 



In this case, the boundary conditions play a crucial role in the integration process. In 
general, considering the inverse relation W^B ^ ^ — A, we have a new set of equations 



where t// := 

2a2 + Oi 
^/7/(x, x'), we have 



2bi(3a2 + 61) ^4^^ ^ 

ai{2a2 + bi) 2a2 

Oi 2a2 



602 + 61 2 



61 



ai 



2a2 + bi 



A = SvrGp - Tii , (A.2) 
(A.3) 



2a2 + bi 

A{^)^Yii y"dVa/Kx,x')^//(x') 



V $. Now by introducing a new Green function 



(A.4) 



where 

a//(x) 87rGp(x) - T7/(x) , (A.5) 

and Yji is a constant having the dimension of a length^ ^. Eqs.([3T|l - ((32|) become 

26i(3a2-|-&i) 4 , , 6a2 + &i 2^ . /n 
ai(2a2 + bi) 2a2 + bi 



ai 



-g,Kx,x') = r2-i5(x-x'), 



(A.6) 

2a2 + Oi 

V2$(x)5,b-$(x),,b = 0. (A.7) 

The general solutions of Eqs. ([?Tl) - (|32p . by introducing the Green function C///(x, x') 
and the function $(x), are 

8a2 + 36i 



C/(x) 



2ai(2a2 

Qi 
2a2 + 61 
1/7 



^17, / dVe„(x,x') 

oi) 7 



87rGp(x') 



*(x') 



2a2 



87r(2a2 + 61 



2a2 + bi 
d^x'd^x 



V^,<i>(x') 



i3„/j3„//5//(x',x") 



2a2 + 61 
1 



$(x") - 



2a2 



|X - X' 
72 



87rGp(x") 



2a2 + 61 



Vi„<i>(x") 



87r(2a2 



, <i>(x') 



(A., 
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T/(x) 



(4a2 +61) 
2ai(2a2 + b 
ai 



-^Yii J d^x'gj/(x,x') 



2a2 + bi 
Yii 



2a2 



47rai(2a2 + 61) 
ai 



2a2 + bi 
1 



<i>(x") 



2a2 + 61 

d^x'd^X' 

2a2 



V2,$(x') 



87rGp(x") 



2a2 + bi 



V^„$(x") 



.3 H^') 



47rai(2a2 + &i) 7 " " |x - x'| ' ^^'^^ 

With the second relation between A and B, Ea. (jA.2[) is a fourth - order equation, 
but in this case, the potentials U, V are linked to A through repeated integrations of 
(jA.Sp - (|A.9P : for the first choice, we have only the integral (I44|) - (I45p . If we consider, 
instead, Eq. (|A.6p . we can find a similar Green function Qii(x,x') for the solutions 
([X8 l) - (|X9|l . In fact, for r 7^ 0, we have 

26i(3a2 + foi)Vf ^//(r) + ai (603 + 6i)V2g„(r) - a?5//(r) = , (A.IO) 

and its solution is similar to (l47l) : 

1 



Gidr) 



Kii ^e V ^1 +Kii^ie\ 
Kii,ee^/^^^^^'' 



KiirC 



2(3a2+i>l) 



(A.ll) 



where, as above, -f^z/.s, Kn^^, Kn ^^ Kh q are constants. It is worth noticing that it 
is not possible to factorize the Laplacian in (|A.6P and, in terms of Fourier transform, 
a vanishing pole is not present. If the Fourier transformed function has no pole in 
the origin (fc^ = 0), this means that, in the field equation, we do not have all terms 
containing a Laplacian but, some term can be interpreted as the mass for the field. On 
the other hand, if it is not possible to factorize a Laplacian in the field equation, this 
means that the pole /c^ = is absent in the Fourier transform of the Green function and 
a Newtonian potential scaling as r~l is not present. Let us remember that a potential 
scaling like 1/r is used in the Green function proportional to / d'^kk"^ e*'' '^"'' ^ 
Furthermore, the analogy between the two approaches is complete when we consider 
the link between them being: 

e/7(x,x') = V2g/(x,x'). (A.12) 



Appendix B. Conventions & Dimensions 

In order to fix the notation, we provide two tables with definitions. Table IBll gives 
an overview of the geometrical quantities used in the paper. We make use of the 
summation convention over identical upper and lower indices throughout the paper. 
The dimensions of the different quantities appearing throughout the work are displayed 
in Table |B2l 
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Object Definition / Convention 



Index ranges a, /3 = 0, 1,2,3; i,j = 1,2,3 

Flat metric r]a0 = diag(l, —1, —1, —1) 
Coordinates = (x'^ , , x"^ , x^) — {ct , x^ , x"^ , x^) 

Vectors v — {v^,v^,v^); V — {d/dx^,d/dx'^,d/dx^) 

Symmetrization T(^a\i3...'',\5) = ^iTaf:)...y5 + Tsfj,,,^a) 

Kronecker (5^ = 1 if a = 0, else 

Connection T^^ = ^g"^" (g^,a.v + 3;/^.^ - 9tiv,a) 
Riemann tensor = r^,_^ - r^^_, + T^.F^^ - T^^F^, 

Ricci tensor = ^av 

Table Bl. List of conventions and definitions. 



Dimensions of Physical Quantities 



gal3, Sa/S, Q, Ol, 



r^iy ^, RaPfiiy ^ , Raff ^ , Yl , Yj/ , , ' 

S-i, bj, Ki^2, Ki^3-\ Ki,4-\ Ki,5-\ 

Ki^7~'^, Kii^3, Kii^i, Kii^5, Kii^Q, r 



Kg M 



kg 

kg_ 



c, c/2, y2, /ii^, u«, ns, A2, $2, r/ 

B 

p 

G 



Table B2. Dimensions of the quantities considered in the paper. 
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